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Abstract 

The energy spectrum of nucleons in high-density nuclear matter 
is investigated in the framework of relativistic meson-nucleon many- 
body theory, employing the 1/N expansion method. The coupling of 
the nucleon with the particle-hole excitations in the medium flattens 
the spectra in the vicinity of the Fermi surface. The effect grows loga- 
rithmically for increasing density and eventually leads to instability of 
the normal state. The validity of the mean-field theory at high density 
is criticized. 
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One universal character shared by most normal Fermi liquids is the en- 
hancement of the effective mass around the Fermi surface 0. Typical indi- 
cations observed empirically are the density of single-particle levels around 
the "Fermi surface" in nuclei, and the logarithmic (T 3 lnT) correction to 
the specific heat at low temperature T in normal liquid 3 He. It is now well 
established that these phenomena due to the energy- dependence of the effec- 
tive mass reflect the dynamics beyond the mean-field theory: The relevant 
mechanism is the coupling of the single-particle motion to other excitation 
modes, e.g., particle- hole (ph) excitations in the case of nuclear matter. As a 
result, the nucleon effective mass at normal density is close to the free mass 
M or slightly larger at the Fermi surface while it is about the traditional 
value (~ 0.7 M) away from the Fermi surface. 

The effective mass at the Fermi surface determines the density of states at 
the Fermi surface. Therefore, the dynamical couplings to the single-particle 
motion could influence the observables significantly. In this light, the density- 
dependence, especially the high-density behavior, of the effective mass would 
be relevant to descriptions of high-density matter, needed to understand, 
e.g., stellar structure or evolution and high-energy heavy ion reactions. 

The investigation of the high-density behavior of the effective mass re- 
quires a relativistic framework of Fermi liquids. In this paper we compute the 
quasiparticle spectra including higher-order many-body correlations employ- 
ing a relativistic many-body theory: In order to exhibit clearly the points 
avoiding complicated aspects such as the compositeness of hadrons and the 
possibilities of new degrees of freedom at high density, we use the simplest 
relativistic model of nuclear matter, the "era; model" ||: 

1111 

C = -0 (ifl - gj/) - M + g a a) ip + - {d^af - ^mla 2 - -u 2 ^ + -m£wj, (1) 

with uj^ = d^ujv — dyijj^. ip, cr and uj^ are the nucleon, a- and omeson fields. 
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Since this model is renormalizable, it allows a systematic investigation of 
the many-body correlations. One convenient scheme for this purpose is the 
1/N expansion where N is the extended number of nucleon species [||. The 
expansion can be obtained by applying the rescaling: g a — > l/y/Ng a , g u — > 
1/yfNgu] <J — »■ y/Na, uo^ y/Nuo^ in the lagrangian ([[]). The energy density 
of nuclear matter up to the next-to-leading order is given by £ = £ + £ e + £ c 
HQ. £q denotes the Hartree term which is the leading order of O(N), while 
£ c and £ c denote the exchange and the RPA-type correlation energy terms, 
both of which are of 0(1). The formulae for these terms can be found in 
ref. 0. Based on this framework, a successful description of the saturation 
property, the equation of state, and the Fermi-liquid properties of nuclear 
matter has been obtained 0, |], 0], improving the Hartree description. The 
essential ingredients are the many-body correlation effects due to £ c . 

The quasiparticle energy corresponding to our energy density can be ob- 
tained by taking the functional derivative with respect to the Fermi distri- 
bution function n(p) following the Landau theory of Fermi liquids |9| and its 



relativistic extension 1C 
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<P) = = eo(p) + ee(p) + e c (p), (2) 

where the three terms derive from the corresponding three terms of our energy 
density. The familiar Hartree term eo, which is of 0(1), contains the scalar 
and the vector mean-field potential contributions. The other contributions 
of 0(1/N) are given by 

M* _ 

e x (p) = -^-^u(p)S x (p)u(p) (x = e, c), (3) 



where E*(p) = \Jp 2 + M* and M* = M - g a a with a the ground-state 
expectation value of the cr-field. u{p) is the Hartree spinor normalized by 
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uu = 1. Z7 X is the nucleoli self-energy shown in fig. |T]; x = e and c correspond 
to the first (exchange) and the following (RPA) graphs, respect ivelyj] 

The counter terms to renormalize the divergences of the self-energies are 
determined through eq. @ consistently with the renormalization conditions 
for the energy density: We renormalize all vertex parts in free space at the 
scales q 2 = and $ = M for the meson and the nucleon legs, respectively. 
The Landau ghost singularity appearing in the one-loop meson propagators 
is removed by the method based on the Kallen-Lehmann representation pro- 
posed by Redmond and extended to finite density in refs. |j, |, [|. 

If the nucleon Hartree propagator connecting the external legs in the 
diagrams of fig. [1] is separated into the "Feynman part" (F) and the "density- 
dependent part" (D) , we obtain the conventional separation of I7 X into Z" x f 
and U x b- For the present purpose, it is convenient to employ a new separation 
of the self-energy based on a Wick rotation (WR) of the loop integration in 



Zxp @, PH : e x (p) = e x w(p) + e x p(p), where e x w involves the Wick rotated 
euclidian 4-integral, and e x p contains an on-shell intermediate nucleon line, 
arising from the nucleon pole due to the WR and the contribution of E x n. 
e x p can be expressed compactly: 

e xP (p) = 4 fj ^ £ dcos 9 Up, q), (4) 

with 9 the angle between p and q. / x (p, q) are the exchange-type interactions 
between two quasiparticles with momenta p, q. f e involves the noninteracting 
meson propagators, while f c takes into account the modifications due to the 
one-loop meson self-energies. 

The slope of the quasiparticle spectra at the Fermi surface defines the 



^or convenience of the following discussion, we employ a slightly different grouping of 
the next-to-leading terms compared to the preceding works |], |j] . 
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Fermi velocity: 



de(p) 

v F 



d\p\ 



+ v Fe + v Fc , (5) 



where the three terms derive from the corresponding terms of eq. (§). In 
the NR treatment, v F = p F /M c g with M c g the effective mass at the Fermi 
surface. If, as we will show below, the enhancement of M e g becomes stronger 
and stronger with increasing density, v F can eventually become a decreasing 
function of the density. In particular, v F — > means that the spectra become 
flat around the Fermi surface. 

In fig. H] we show v F as a function of p F . The parameters are taken from 
refs. |3|, |J and are shown in table [TJ: The dotted line shows the Hartree 
approximation, i.e., the first term of eq. (pj) using the "Hartree" parameter 
set; it approaches the causal limit "1" at high density. In contrast, the solid 
line, which shows the full result of (IS) using the parameter set B, crosses 
v f = at the critical p F — pp, and v F becomes negative for higher densities. 
The dashed line, including only one density-dependent bubble instead of 
the full RPA propagator^] shows a behavior similar to the full case, but p F 
becomes smaller. We also show the full result of (0) using the parameter set 
A by the dot-dashed line; the comparison with the solid line demonstrates 
the sensitivity of p F to the parameters (see eq. (|ll]) below). In fig. |] we show 
the quasiparticle energy of eq. (0) for the case of the solid line of fig. |2|. For 
high densities, the spectra show the "anomaly" around the Fermi surface, i.e., 
they are flat or decreasing due to the inclusion of the higher-order self-energy 
terms. f\ Thus, the effect grows with increasing density. 
2 The Feynman bubbles are summed up in this result. 

3 The latter case with the negative does not correspond to the normal state and there- 
fore, strictly speaking, the Landau theory is not applicable. This point will be discussed 
later. 
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From fig. 0(b), we see that this effect is due to e c p, which is a decreasing 
function around the Fermi surface. (This is seen already at normal density, 
but for higher densities this effect becomes more pronounced and shows up 
in the total results.) The decrease of e c p around the Fermi surface could be 
expected from eq. (§]) if the result of the angular average of f x (p, q) assumes 
a positive value (see also ref. [IT] for the NR case). 

Now can we understand the origin of the growth of the anomaly due to 
e c p? For this purpose it will be useful to examine the leading high-density 
behavior of vp. We set for pp — > oo vp — > 1 + 5vp e + 5vp c = 1 + 5vp, where 
"1" is due to the first (Hartree) term of eq. (|5|) while 5vp x (x = e, c) denote 
the asymptotic limits of the following two (0(1/ N)) terms Vp e , vp c : 



Svp x » -4 

7T J— 



£dcos0 Up,q)\ p=q=pF . (6) 



The r.h.s. is the result obtained by substituting e x p of eq. (§) into eq. 
"w" means that the contribution due to e x \y is neglected because of the above 
discussion of fig. |||. (In fact, it can be shown analytically that e x w does not 
alter our conclusions JTljfl.) 

It is convenient to separate the integrand of eq. (Q) according to the types 
of the possible meson modes in the medium |4j]: 

UP,q)\ P=q=PF = -^EA^cos ^(cosfl), (7) 

opp j 

where I = a, M, L and T, corresponding to the a-, mixed, longitudinal and 
transverse contributions. h\ is a dimensionless function accounting for the 
scalar and the vector vertex structures. A e \ are the exchange contributions 
due to the noninteracting a- and cu-meson propagators while A c \ give the 
modifications of A e i due to the RPA correlations. 

To extract the leading asymptotic behavior of (0), we set M* = 0; because 
there is no singularity for M* = and from the dimension counting, any 



6 



contributions omitted here will be suppressed for p F — > oo compared to the 
retained ones. In this limit it can be easily seen that the contributions due 
to the M-mode are suppressed compared to the others. 

In this case, the relevant contributions (I = a, L and T) have the form: 



9i 1 (Mt) 



(8) 



2t + K( \t + K{J 1 

where t = 1 — cos 9, g L = g T = g^, m L = m T = m u , nj = mj/2p F , 
and fli{t) = IJi{p — q)\ p=q=PF /2p| with III the 1-loop self-energies for those 
modes.0 a\ give the contributions involving n bubbles for the I-mode, such 
that A el = a? and A cl = Y, n >i(- l ) n ~ la i- 

If the integrand of eq. (j^) given by eqs. (|7J), were finite everywhere 
in the integration region, the result would coincide with the naive dimension 
counting which leads to 5v Fx — > const as p F — > oo. For p F — > oo, however, 
k\ — > and therefore the integral diverges at = 0. This divergence implies 
that 5u Fx could contain positive powers or logarithms of p F . (It can be shown 
13 1 that there is no other singularity.) — > corresponds formally to zero 



meson masses; therefore, these singularities for # ~ 0, which now prove to 
give the most dominant contributions for p F —>■ oo, correspond to the infrared 
singularity in the massless limit. 

By substituting the analytic formulae of hi @] and ili [§, into eqs. 
(@)-(@) and by extracting the singularities for 9 ~ 0, we can obtain Sv Fx up 
to logarithmic accuracy. For x = e, we obtain 

= (gl x 0(1)) + ( 1 -(^) 2 U4 + 0(1))) 



2n J [ rri^ 

-i m 2 l^+o(i)\). " (9) 



2 V2tt/ m 2 , 



4 In this work 77l denotes the minus of "27l" m refs. [|[ ||. 
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The three terms show the contributions due to the relevant three modes a, 
L and T, respectively. Though both the contributions due to the L- and 
T-modes grow logarithmically, reflecting the logarithmic infrared divergence 
at 9 — > 0, they cancel out, leading to 6v Fe = const. The contribution to Sv Fc 
due to one bubble insertion, i.e., due to the n = 1 term of eq. (H), can be 
obtained similarly: 

5v Flh = fexO(ln^)) + f-f^Vl^-^n||ln4 + 0(lng)ll 



2ttJ I m 2 , 3 M 2 m 2 

\ \ CO to s / 

The coupling to the ph excitations makes the infrared divergences stronger;^ 
the most dominant contribution is due to the linearly divergent integral due 
to the L-mode, and the result grows like p|. However, it can be easily 
seen that the insertion of more bubbles causes stronger infrared divergences, 
leading to terms of higher powers in pp. 



This situation is reminiscent of the high-density NR electron gas [13 



This suggests that we have to sum up all the bubble graphs, forming the 
complete RPA series. Therefore, the infrared structure of vp due to the cou- 
pling of the ph excitations naturally leads to our 1/N expansion scheme. By 
summing up all the bubbles, we obtain for the total Sv^: 

Sv F = (g a ( PF ) 2 x 0(1)) + (Upf) 2 x 0(1)) 

The resummation of the ring graphs renders the infrared divergence mild, 
leaving only the logarithmically growing term due to the T-mode (compare 
the dashed and the solid curves of fig. H). Thus the decrease of v F , which led 



5"V 



h" includes the negative-energy states. 
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to the inversion of the spectrum of fig. |3], is due to the logarithmic correction 
flllD due to the self-energies beyond the mean-field theory. 

The resummation also replaces the original coupling constants gi renor- 
malized at the meson momenta q 2 = by the running coupling constants 
<7i(Pf) a t the scale pp [fL5|l : Before applying Redmond's method to avoid the 
Landau ghost H, they are related by the one- loop formulae: for pp M, 
9u;(pf) 2 ^ gl /7l - |§r In ff^ , and similarly for g a (p F ). £i(pf) of eq. ([11]) 
should be understood as the result of Redmond's method; they grow for in- 
creasing pf toward finite values for the bare coupling constants Thus the 
growth of the last term of eq. (|11]) is enhanced by the growth of ^(pf)- 

Our result can be summarized by the behavior of the screening mass in 
the medium. The infrared behavior of the I-mode contribution is governed 
by the screening mass (squared) m 2 cl , which is obtained from the self-energy 
i7i(g) by letting q° -> first followed by q -> (see eqs. ©-©). The 
Lorentz invariance guarantees m 2 ch = m 2 cT , which would lead to the cancel- 
lation between the L- and T-mode contributions (see eq. @). However, the 
coupling of the ph excitations breaks the invariance: The L-mode acquires a 
nonzero Debye screening mass, which is m 2 ch = 2g 2 pp/ir 2 for p-p — > oo, but 
m scT = d ue to baryon current conservation [0. As a result, the cancellation 
between the two modes becomes incomplete, leaving the result of eq. (|TT|). 

Now we discuss physical interpretations of our results. 

Above the critical p F where vp = 0, the normal state is unstable. In our 
calculation of the ground-state energy density and the quasiparticle energy 
we assumed the normal state 0. If Vp < 0, however, this assumption is not 
valid. This suggests a stability condition: v p > for a normal Fermi liquid, 
in addition to the wellknown conditions for the dimensionless Landau- Migdal 
parameters 0. If the ground state is not normal, one should re-compute the 
energy density self-consistent ly allowing the possibility of new configurations 
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like the so-called "Fermi-gap" state: If one literally accepts the results of fig. 
[3], this indicates a 1st order phase transition at p-p = pp from the normal 
state to the Fermi-gap state.0 In contrast to the Fermi sphere for the normal 
state, the nucleons occupy the momentum states within the inner sphere and 
the outer shell for the Fermi-gap state. In our model, pp = 4.07fm _1 which 
would imply the critical density for neutron matter g c ~ 15^0 • 

We also note another possible scenario: Near the critical point vp = 0, the 
terms of higher order than we considered here, which would give contributions 
~ (^(Pf) 2 lnp|) n {n > 2), also could be important. It would be a very hard 
job to compute these contributions. Intuitively, one could regard the result of 
eq. ( |TTD as the leading approximation of vp ~ 1/ (l + (^(pf)/27t) 2 /2 x In (pp/m^) 
In this case v f = would not occur for any finite pp; vf — > asymptotically. 

In connection with these scenarios, we stress that the (formal) high- 
density limit based on the lagrangian (|l|) does not correspond to the mean- 
field theory, contrary to the usual claim 0: For the first scenario this is 
obvious; the second one implies an infinite density of states at the Fermi 
surface for pp — - > 00, which never be reached in the mean-field theory. The 
usual claim tacitly assumes the absence of infrared singularities, which reflect 
large long- wavelength fluctuations of the fields, from the higher order correc- 
tions. This assumption leads to the validity of the naive dimension counting, 
and therefore to the dominance of the mean-field contributions because they 
contain the largest number of independent fermion loops in each order of the 
interaction. As discussed above, the effects which lead to the modification of 
this naive argument are due to the strong infrared behavior induced in the 
high-density medium. 

6 The possibility of the Fermi-gap state for high-density nuclear matter was pointed 
out in ref. [Q in the Hartree-Fock approximation of the lagrangian (|l|). However, as 
stressed above, the mechanism in our case is beyond the Hartree-Fock approximation 
(recall Svf c = const). 
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This conclusion for the high-density behavior seems to be inevitable for 
relativistic descriptions: The relevant infrared behavior is due to the ph 
excitation in the vector channel (see eqs. (P)- (|TT|) ).p1 However, the vector 
boson degrees of freedom are indispensable to keep matter from collapsing 
at high density by supplying the repulsion. We also note that the conclusion 
might be valid even if the compositeness of hadrons were taken into account: 
Even though the vertex form factors were introduced, they could not influence 
the results dominated by the contribution from the infrared region (q 2 = 0). 

Though we have shown the possibility of the Fermi-gap state, the tran- 
sition density g c appears to be above the critical densities to more familiar 
new phases like quark matter phase, pion condensed phase, etc. However, 
our result was obtained just within a simple model, and g c is very sensitive 
even to small corrections, which might be due to higher orders in 1/N or 
other hadronic degrees of freedom like ir, p, etc. After the inclusion of these 
effects, the competition or coexistence of the Fermi-gap state with more fa- 
miliar new phases would be an interesting aspect of high-density hadronic 
matter. 

In conclusion, we have investigated dynamical effects in high-density nu- 
clear matter due to the higher order many-body correlations. The anomaly 
in the fermion spectra, corresponding to the famous effective mass enhance- 
ment in NR Fermi liquids, grows with increasing density. The effect is due to 
the strong infrared behavior caused by the coupling of the ph excitation to 
the single-particle motion. The picture of high-density matter emerged from 
our investigation is much more complicated than the mean-field description. 
Though a more realistic treatment would be required to fix the relevance to 

7 In view of this, the extension of our approach to gauge theories would be appealing, 
and is under investigation for cases of cold QED and QCD plasma with finite fermion 
density fl5j . 
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phenomenology, the universal character of this phenomenon already poses 
various interesting aspects for high-density matter, like the possibility of the 
Fermi-gap state, the thermodynamic properties and transport phenomena in 
this phase, and the extension to gauge theories. These points will be reported 
in future publications |15[ . 
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Figure captions 



1. Feynman graphs of the nucleon self-energy of 0(1/N). The dashed lines 
in the figure denote the noninteracting a- and a>-meson propagators; the 
solid lines denote the nucleon Hartree propagators on the background 
meson fields and the Fermi sea. 

2. The Fermi velocity as a function of the Fermi momentum. For expla- 
nation of the curves, see text. 

3. The real parts of the quasiparticle energy (eq. (0)) are shown as func- 
tions of momentum by the solid lines, (a), (b) and (c) are for cases 
of pf = 1.30, 4.07 and 4.50fm _1 , which respectively corresponds to 
vp > 0, = and < 0. The dot-dashed line shows the Hartree con- 
tribution contained in the full result. The insertion in (b) shows the 
0(l/iV)-contributions as functions of momentum. The dot-dashed line 
shows e e of eq. (|2|), while the solid line shows e c . The dotted and 
dashed lines are the separated contributions e c w and e c p contained in 
the solid line. 

Table captions 

1. The parameters used for numerical computation. 
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Table 1 





g 2 j4ir g 2 j4ir m a [MeV] m w [MeV] 


Hartree 

l/N A 
B 


6.23 8.18 550 783 

2.24 3.65 550 783 
2.60 4.22 650 783 
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